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Abstract 

After a pedagogical introduction to the calculation of weak matrix elements on the 
lattice, I review some of the lattice's most recent predictions concerning the weak decays 
of i?-mesons. Amongst the topics covered are the determinations of the leptonic decay 
constant /b, of the form factors relevant for semi- leptonic B D{D*)ii> decays and 
of the hadronic matrix element which describes the rare decay B — > K*^. Emphasis 
is placed on the results of the UKQCD Collaboration. Corrections to the heavy-quark 
limit are discussed extensively. 
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1 Introduction 

Weak decays of hadrons are a very rich source of information about the Standard Model. Not 
only do they enable us to determine from experiment the parameters which are associated 
with the flavor sector of the Standard Model but they also provide a rich testing ground for 
understanding the non-perturbative dynamics of the strong interaction. These two aspects 
of hadronic, weak decays are in fact inseparable. To be more specific, the Standard Model 
prediction for the rate of such a decay follows the pattern 



from which it is clear that in order to extract the Cabibbo-Kobayashi-Maskawa (CKM) factor 
from an experimental measurement of the rate, one has to understand or at least calculate 
the non-perturbative QCD factor. The lattice formulation of QCD, together with large 
scale numerical simulations, provides a natural framework for evaluating this factor. It is in 
fact the only systematic, first-principle approach we know for quantifying non-perturbative, 
strong-interaction dynamics with the QCD lagrangian as a starting point. This should not 
be taken to mean that it is only approach or even the most efficient approach, as QCD 
sum-rules, low energy effective theories and quark models have provided many important 
and reliable results throughout the years. 

The latest chapter in the story of hadronic weak decays concerns hadrons containing a b- 
quark. The decays of these hadrons are interesting phenomenologically because their study 
will enable us to determine the least well know column of the CKM matrix. They are also 
interesting theoretically because they provide a ground for testing heavy-quark symmetry. 
Heavy-quark symmetry is a symmetry of QCD that arises in the limit that the mass of the 
heavy quark is much larger than the QCD scale Aqcd- In that limit one finds that the 
dynamics of the light quarks and gluons coupled to the heavy quark become independent of 
the heavy quark's fiavor and spin. To the extent, then, that Aqcd is negligible compared 
to the masses of the charm and beauty quarks, QCD exhibits a new SU (4) spinx flavor sym- 
metry which acts on the multiplet (c t,c |,6 t,6 |)[||, ||. This symmetry tremendously 
simplifies the description of the decays of hadrons containing a heavy quark. It has in fact 
been incorporated into the framework of an effective theory known as Heavy Quark Ef- 
fective Theory (HQET). This makes it possible to systematically calculate corrections to 
the symmetry limit order by order in inverse powers of the heavy-quark masses much in the 
same way that Chiral Lagrangians enables one to calculate corrections to the predictions 
of the spontaneously broken chiral symmetry of QCD in the light quark sector.0 I would 
like to emphasize, at this point, that the lattice is supremely well suited for studying the 
range of applicability of heavy-quark symmetry, for the masses of the heavy quarks are free 
parameters in lattice calculations and the dependence of results on heavy-quark mass can 
be studied in detail. 

"'^For a comprehensive review of HQET and heavy-quark symmetry see the lectures of K. Zalewski in this 
volume or the review of 0. 
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In the present article I will review some of the latest results of lattice calculations of weak 
matrix elements of beautiful hadrons. I will concentrate mainly on the results of the UKQCD 
Collaboration because this article is not meant to be thorough and systematic review of the 
subject but rather an introduction to the methods, possibilities and limitations of this rapidly 
growing field. In Section ^ I will provide a introduction to lattice calculations of weak matrix 
elements. For a more complete review, see the lectures of R. Gupta in the present volume 
or the reviews in [Q. For reviews about 6-physics on the lattice, see Ref.Q. In Section 
^ I will present results for the 5- meson decay constant, fs, and will describe what these 
results teach us about heavy-quark symmetry. I will then turn to the subject of semi-leptonic 
and — i> Z)*g) decays in Section |[ I will present determinations of the Isgur- 
Wise functions C,u,d and C,s which parametrize the strong interaction effects to leading order 
in heavy-quark mass in these decays. I will show how the function C,u,d can be used to extract 
the CKM matrix-element Vcb from an experimental measurement of the differential decay 
rate for B — > D*ii> decays. I will also present many tests of the heavy-quark symmetry 
which involve comparing B ^ D and B ^ D* decays for many different values of the b and 
c quark masses. Section || will be dedicated to a lattice evaluation of the rate for the rare 
process B —>■ K*'~f. This process occurs through a flavor-changing-neutral-current, penguin 
diagram and was first observed in 1993 by the CLEO collaboration!^. I will end the article 
in Section |^ with a summary and conclusions. 



2 Matrix Element Calculations on the Lattice 
2.1 What We Calculate 

To obtain matrix elements in lattice calculations we compute expectation values of products 
of gluon and quark fields at different spacetime points. We do so using the path integral 
formulation of QCD in euclidean spacetime. Thus, given a product 0{xi, ■ ■ ■ of quark 
and gluon fields, we compute n-point functions: 

• • • , x„)) = ^J dA^d^pdij e-^«^^0(xi, ■ ■ ■ , x„) , (2) 

where Sqcd is the QCD action and 

Z = J dA^dtpd^j Q-^QCD _ (^3) 

The simplest example of an n-point function is the 2-point function. To determine the decay 
constant of a 5-meson, for instance, we would evaluate the 2-point function 

^(g7VKx,t)&7'g(0)) , (4) 

X 

where g is a light-quark field and the sum over x ensures that the meson is at rest. If we insert 
a complete set of states between the two operators in Eq. we find that this correlator 
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decays exponentially in time in the limit of large t. The rate of the decay is governed by 
the mass of ground state and its amplitude is proportional to the weak matrix element, 
(0|g7'^7^6(0)|S), we are after: 

The next level of complexity is the 3-point function. The matrix element {D{p')\c'~f'^b{0)\B{p)) , 
for instance, is obtained from the 3-point function 

E e-*'-^-^'^-"(c7'g(y, tf) c^K^, t) b^qm (6) 

where q = p — p'- In the limit that t, tf — t —>■ oo, this 3-point function reduces to 

(0|g7^6(0)|D(pO) (D(pO|c7^&(0)|g(p)) (i?(p) 1^7^^(0)10) ^_^^,_^^(,^._,) ^ 

4Ee)Eb 

where Eb and Ed are the energies of the B and D mesons, respectively. 



2.2 How We Calculate 

The first step in evaluating the path integral of Eq. (^) is to approximate spacetime by a 
finite, hypercubic lattice. Quark fields are then placed on the sites of this lattice and the 
gluon fields on the links between these sites. This reduces the integral over gluon and quark 
fields to an integral over a finite number of degrees of freedom which makes it amenable to 
numerical methods. It is important to note, at this point, that to preserve gauge invariance 
on a discretized spacetime one works with elements of the group-the link variables U^{x) = 
e*''^"^'"^^^ -instead of with the gauge fields A^{x) (see Ref. for example), i This means, for 
instance, that the integrals over gauge fields in Eqs. (H) and (0) are actually integrals over 
the compact link variables f/^(x). 

As a first step in the evaluation of this lattice path integral, we generate a number of 
gluon configurations by a Monte-Carlo simulation[|] with the probability distribution 

e-'^s'""" X det(^-h M) 
V ^ ^ , (8) 

where Sgiuon is a discretized gluon action|^, a discretized SU{3) covariant derivative 
in the fundamental representation and M, the quark-mass matrix. Then, for example, to 
obtain the propagator of a 5-meson, 

{u-f%{x)b-r^u{0)) , (9) 



^Here g is the gauge coupling constant and a, the lattice spacing. This lattice spacing acts as a non- 
perturbative regulator: modes with momenta larger than n/a are absent regardless of whether the theory is 
used perturbatively or non-perturbatively. 
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we use Wick's theorem to reduce the quark bihnears in Eq. (||) to a product of quark 
propagators and do so for each one of the gluon backgrounds. The meson propagator 
then becomes the following average of a product of quark propagators and Dirac matrices 
over gluon backgrounds, 

{u^'b{x) h'um = ^ E Tr {I'S.iO, X- U.h'S^ix, 0; U,)) , (10) 

which can be evaluated numerically once the propagators, 5'b(u)(0, x; f/^), for the b and u 
quarks in the gluon background f/„ have been calculated using matrix inversion algorithms 



such as the conjugate gradient algorithm |110[ . 



2.3 Limitations and Sources of Errors 

Even though the numerical evaluation of the path integral of Eq. described in Section 



0| is an ab initio calculation starting with the QCD lagrangian, the final results have 
errors which are due to the approximations that have to be made. These errors can be 
classified into two broad categories. There are statistical errors which arise because we use 
statistical methods to evaluate the path integral and systematic errors which arise because we 
approximate spacetime by a finite lattice of points and because we usually have to neglect the 
contributions of fermion loops. One may wonder, at this point, what advantage the lattice 
has over other means of evaluating hadronic properties. The difference lies in that the errors 
made in a lattice calculation can be reduced systematically by increasing the number of 
configurations and the physical volume of the lattice, by decreasing the lattice spacing or by 
designing lattice actions which converge to the continuum limit more rapidly. 



2.3.1 Statistical Errors 

Statistical errors arise because we approximate the integral over gauge fields, which is an 
integral over an infinite number of configurations, by a sum over a finite number, N, of 
configurations. According to the central limit theorem, the error made on n-point functions 
is proportional to in the limit of large A^, so that statistical errors decrease as the 

number of configurations increases. Typical values of A^ in QCD calculations are on the order 
of 100 and the corresponding errors are on the order of a few percent for many quantities. 
We evaluate these errors using standard statistical methods such as the bootstrap or jacknife 



method |ri|]. These methods enable us to approximately determine the distribution of results 
we would find were we to repeat our A^-configuration calculation many times. The statistical 
error we quote on our result is chosen to include to the central 68% of this distribution. 
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2.3.2 Discretization Errors 

Discretization errors are due to the fact that we approximate spacetime by a discrete lattice 
of points. To see how these errors arise, consider the symmetric lattice derivative of a 
function, /, of spacetime in the limit that the lattice spacing a is taken to zero: 



where /i is a unit vector in the /x-direction. It is clear from Eq. (|TT]) that approximating 
the derivative of / by a finite difference introduces an error of order a? . Before exploring 
these discretization errors in more detail, though, I shall first describe how we determine the 
lattice spacing a from our simulations. 

Setting the Scale 

As you well know, the bare coupling constant g{a) is related to the cutoff a through dimen- 
sional transmutation. Before performing our simulation, however, we only know what this 
relation is through perturbation theory up to some finite order in g{a). So it is preferable 
not to fix both Quantities from the start. What we do is pick a value for the bare coupling, 
(5 = 2Nc/g'^{a), B and rewrite the QCD lagrangian in terms of dimensionless quantities, mea- 
sured in units of the lattice spacing a. All explicit dependence on a then disappears from 
the problem and a need not be fixed before performing the simulation. Once the simulation 
has been performed, however, a must be determined so that physical dimensions may be 
restored. To do this, we pick a dimensionful, physical quantity such as the string tension a. 
The lattice spacing is then the ratio of the experimental value, ct'^^p*, to the dimensionless 
quantity, o"^^**', that our simulation yields, raised to the appropriate power: 



To set the scale any dimensionful quantity will do though different quantities will lead to 
more or less accurate determinations of the lattice spacing, i In Table |l|, I show the values for 
that our collaboration obtains from various physical quantities. All of these quantities 
should, in principle, yield the same lattice spacing assuming, of course, that QCD is correct. 
The reason why the values in Table |I| do not agree within the quoted statistical errors is 
because we also make systematic errors. In fact, the range of values that we obtain for 
is a good indication of the size of the systematic errors we make. 

Having established how one determines the lattice spacing, we return to the discussion of 
discretization errors and how to reduce them. The reason why these errors are particularly 

■^Here, perturbation theory is a useful guide for it tells us approximately what value of /3 yields the desired 
lattice spacing. 

^For quantities that depend on quark masses, one has to first determine the quark's bare mass. 




(11) 




(12) 
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Physical Quantity 


a-'^ (GeV) 


^Ja (string tension) 
vrip (mass of the p) 
rriN (mass of the nucleon) 
rriA (mass of the A) 


2.73(5) 
2.7 t I 
3.0 tl 
2.5 tl 



Table 1: Values of the inverse lattice spacing, a~^, for UKQCD's 
(3 = 6.2, 24^ X 48 lattice as obtained from different physical quanti- 
ties |]12[. Ratios of these different values indicate how well the corre- 
sponding ratio of physical quantities is predicted by the simulation. 
a~^{mp)/a^^{mN), for instance, indicates that the ratio rriN/mp comes 
out 10% low in our simulation but with errors large enough to make 
it consistent with experiment. 



important in simulations of hadrons containing a heavy quark is because these quarks have 
very short compton wavelengths: an accurate description of their quantum propagation 
appears to require a very fine lattice. There are in fact many different ways to describe 
heavy quarks on a lattice. These different approaches fall into two broad categories. In the 
first category, the quark action used is simply a discretization of the full continuum quark 
action. In the second, that of effective theories, this discretized action is further expanded 
in inverse powers of the heavy-quark mass as well as in powers of the heavy quark's velocity 
if the heavy quark is non-relativistic. 

The Wilson Action 

One of the standard discretizations of the continuum quark action is the Wilson actionQ: 

Up{x)q{x + afi) + q{x + ajl){r + 7^)?7^(x)g(a;) 

+g(x)(m„ + ^)g(x)| , (13) 

where r is a number (typically taken to be 1) known as the Wilson parameter. In the naive 
continuum limit, this action reduces to 

Jq(ip+m^)q-^JqD'q + 0{a') . (14) 

Because of the O (a) term, we expect the leading discretization errors, in a simulation of a 
hadron containing a heavy quark, to be on the order of arnq, where mg is the mass of the 
heavy quark. For these errors to remain under control, we must require amq -C 1. Since 
typical values of the inverse lattice spacing are on the order of 3 GeV, the b quark, which has 
a mass of about 4.5 GeV, cannot be simulated directly. Thus, the strategy used to obtain 



S 
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properties of 6-hadrons with the Wilson action is to calculate these properties for a collection 
hadrons whose heavy quarks have a mass in the range of the charm quark mass and then to 
extrapolate these properties in heavy-quark mass to rrif,, using HQET as a guide. But even 
for the charm quark, we expect large discretization errors as arric ~ 0.4 for a typical value 
of the lattice spacing. So if the Wilson action is to be used, one has to find a way to reduce 
these discretization errors. 

The brute force approach to reducing discretization errors is to increase the number of lattice 
sites and reduce the lattice spacing. This requires ever faster computers. The second type 
of approach is to design lattice theories which converge faster to the continuum limit. It is 
to this second approach we now turn. 

Improving the Wilson Action: Normalization of Quark Fields 

The impetus for the improvement that I am going to describe comes from an analysis of the 
free-fermion propagator on the lattice. One finds, using the action given in Eq. (p!3D , that 
this propagator is given, at zero momentum, by 



(15) 



2 

where the term in brackets is the continuum, free propagator for a fermion of mass 

am = In {1 + arrio) . (16) 

In the limit that arrio vanishes, the lattice propagator in Eq. ([15D reduces to the continuum 
propagator for a fermion of mass rrio, as it should. However, when anio ~ 1, as it is for the 
charm, this lattice propagator starts deviating from continuum behaviour. To remedy this 
problem, it has been suggested [P^, Q that one should let m, defined in Eq. (|TBp, be the 



bare quark mass and that one should rescale the lattice quark field q according to 

For the charm quark, this corresponds to rescaling the quark field by a factor of about 1.2, 
quite a large factor. 

Even though this procedure clearly improves the free lattice propagator it must be stressed 
that it is not a systematic procedure and that it is unclear, at this point, what the remaining 
discretizations are once interactions are switched on. For instance, one would expect that 
a^amo and oAqcd discretization errors are still present. 



It has even been suggested recently [ 14 1 that it is possible to use the Wilson action for any 



value of amo-even anio ^ 1-as long as one performs the rescaling of Eq. (p!7D and one defines 



the mass, M, of the simulated hadron as 
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Here again, however, one is entitled to suspect that large a^amo-errors might contaminate 
results for matrix elements. 



Improving the Wilson Action: Symanzik's Procedure 

A more systematic approach to reducing discretization errors was initiated by Symanzik in 
the early eighties [|I^]. Here the idea is to formally remove discretization errors order by order 
in a by introducing higher-order operator corrections to the action. The first step in this 
program, as it applies to lattice QCD, was carried out by Sheikholeslami and Wohlert who 
proposed the O (a)-improved action |I6[ 

gsw ^s"^ - irg^a^ 9 ^m-^a*- ^i^) , (19) 

where P^jy is a discretization of the field-strength tensor F^j,. This action is known as the 
Sheikholeslami- Wohlert action or clover action. The name "clover" comes from the fact that 
a possible expression for P^,^ is the four-leafed-clover product of the four plaquettes that 
originate at x and lie in the //i/-plane. 

It has been shown [P^ ] that when this action is used in conjunction with the field rotation 

g g« = (^1 - y ^) q (20) 

the discretization errors made on matrix elements are reduced from O (a) to O (aga), where 
as is the strong coupling constant. Thus, the discretization errors which we evaluated to be 
on the order of am^ ~ 40% when simulating the charm with the Wilson action on a typical 
present day lattice (i.e. /? = 6.2) are reduced to about a^amc ~ 4% when the clover action 
is used i. 

Improving Lattice Actions: The Perfect Action 



A little over a year ago P. Hasenfratz and F. Niedermayer raised the tantilizing prospect 
of a "perfect action" for lattice QCD. Such an action would yield cut-off independent phys- 
ical predictions on relatively coarse-grained lattices. What they have done is to develop a 
workable procedure for following the renormalization trajectory of an asymptotically free 
theory from the continuum fixed point out to relatively large values of the coupling. They 
have convincingly shown that their procedure produces a "perfect action" for the d = 2, 0{3) 
non-linear cr-model whose range of interactions is short and whose structure is relatively sim- 
ple. They are currently applying their ideas to Yang-Mills theory in four dimensions and are 
finding encouraging preliminary results|T^. It must be said, however, that a perfect action 
can only be used to answer questions which refer to modes which have not been integrated 
out so that their usefulness for simulating heavy quarks may be more limited than it is for 
light quarks and gluons. 



^The value of as used here is the "boosted" value which incorporates the effects of tadpole diagrams 
and is therefore more physical. For the UKQCD simulation at /3 = 6.2 it is as — 3/(27r/3M^), where Uo is a 
measure of the expectation value of a link variable which we take to be Uo = 1/ (SKcrit)- i^crit = 0.14315 is 
the critical value of the hopping parameter. 
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Effective Lattice Field Theories: Static Quarks 

Up to now we have been considering different ways of improving the Wilson action. Another 
approach is to give up trying to treat heavy quarks as fully relativistic fermions and to expand 
physical quantities in inverse powers of the mass, rriQ, of the heavy quark, thereby getting 
rid of large O {amo) discretization errors. The zeroth order in this expansion constists in 



treating heavy quarks as static, spin- 1/2 color sources pO[. The corresponding action is a 
discretization of 

gstat ^ f Q^J^^Q ^ ^21) 



where Q is a two-component spinor field. The propagator for such a quark is trivial as it 
is essentially a product of link variables in the time direction. This makes static quarks 
very easy to simulate, at least in principle. One can show|T^, ^|, however, that the ratio 



of signal to noise for the propagator of a hadron containing a static and a light quark falls 
exponentially and very rapidly in time so that correlators made from static propagators are 
usually extremely noisy. Furthermore, because heavy quarks are static in this approach, 
we cannot study transitions between heavy quarks of different momentum as one must to 
describe decays such as B ^ D*lv. 

Effective Lattice Field Theories: HQET 

In this approach, one allows the heavy quark to have a non-vanishing velocity v and be 
slightly off-shell. To leading order in heavy-quark mass, the relevant action is a discretized 
version of 

gHQET ^ J Q\^^D,Q , (22) 

where Vo = \/TT~v^ and Q is the same two component spinor field as above. Obviously, 
when V = 0, we recover S****"*. 

Like the static action, this action is much easier to simulate than the full action, because 
calulating the quark propagators is an initial value problem. The HQET moreover permits 
the study of transitions between heavy quarks of different momentum which is necessary to 
obtain the form factors relevant for weak b ^ c decays in the heavy-quark limit. However, 
the low signal to noise ratio problems encountered with static quarks are present here too. 

One may also wonder whether this theory is even defined in Euclidean spacetime as the free, 
momentum space propagator has a pole at Po = iv ■ p/vo which may be on the positive or 
negative imaginary axis depending on the relative orientation of v and p. This means that 
the free theory has solutions which grow exponentially in time. One can show, however, that 
the resulting divergences are regulated when these propagators are combined with relativis- 
tic, light-quark propagators in correlation functions for hadrons containing heavy and light 
quarks [p2| . 



Effective Lattice Field Theories: Non- Relativistic QCD (NRQCD) 

In this approach, in addition to expanding the QCD action in inverse powers of the heavy- 
quark mass, one further expands it in powers of the heavy-quark's velocity squared, v^. 
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which is assumed to be small. This double expansion leads, at order Aqcd/"^q and v^, to 
the following action p3|, ^ 



5™^^= /Qt|z)^ + ^ + ^^.B|g, (23) 
J { 2mQ 2mQ J 

where B is the chromomagnetic field and Q is, again, a two component spinor field. 

This action permits one to simulate the b quark directly. It is also simpler to simulate than the 
full theory, for the same reasons S^'^^'^ is. Furthermore, the additional kinetic term in the 
action guarantees that the ratio of signal to noise in correlators which include non-relativistic 



propagators is much better than in correlators which include static propagators pi] . 

One has to remember, however, the NRQCD is an effective theory and that the cutoff 
must be kept finite. Since the expansion which leads to the action of Eq. (p3D improves 
as amq becomes larger, should not be chosen much bigger than about 2 GeV when 
the 6-quark is being considered and should be chosen much smaller when the c-quark is 
involved. This means that one has to begin worrying about discretization errors. Since 
these errors cannot be removed by going to the continuum limit, they must be removed by 
higher dimension operators. Whether such a program can successfully be carried out for 
both spectral quantities and matrix elements has yet to be shown. Moreover, from a more 
practical point of view, the amount of phase space that can be explored in B —>■ D decays, 
for instance, will be severely limited by the fact that one has to limit the size of the momenta 
given to the particles to keep momentum-dependent discretization errors under control and 
for the non-relativistic approximation to remain valid. 



2.3.3 Finite Volume Errors 

These errors are due to the fact that we approximate spacetime by a box of finite volume. The 
first requirement on the size of this box is, of course, that it be large enough to accommodate 
the hadron we wish to study. If this is not the case, volume corrections which fall off with 



inverse powers of the lattice's volume appear[25, 26|. For mesons composed of a heavy quark 



and a light antiquark, the requirement that the box be larger than the meson means that 
aL 3> 1/Aqcd, where L is the number of lattice sites in the space directions, since the typical 
size of these mesons is on the order of 1/Aqcd- 

The second requirement, in the case of peridic boundary conditions, is that the box be 
larger than the range of the strong interaction so that the hadron cannot interact with its 
many copies[^, Since this range is determined by the mass of the lightest hadron-the 
long range component of the strong interaction is governed by pion exchange and falls off 
exponentially with distance as e""^'"''- we must ensure that copies are at the very least one 
pion compton wavelength away. Because there is a limit on the size of the lattices that we can 
work with, this requirement is usually turned around and used to constrain the mass of the 
light quarks. If we want our finite volume errors to be smaller than about 5%, neighboring 
copies must be at least three pion wavelenghts away (e~^ ~ 5%). This means that the mass 
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of the lightest allowed quark must be such that a pseudoscalar meson which contains this 
quark and its antiquark has a mass mp with 



3 . 1 

<aL- 

mp Aqcd 



(24) 



For the lattice used by the UKQCD collaboration {a~^ ~ 2.7 GeV and L = 24) we would 
therefore require that the lightest pseudoscalar meson have a mass larger than ~ 600 MeV 
(Aqcd — 250 MeV). This constraint can, in fact, be slightly relaxed since we work in the 
quenched approximation.!^ 

Before leaving the subject, I would like to point out that our collaboration has recently 
uncovered an unexpected finite volume volume effect [^] in trying to implement a new idea 
for obtaining the slope of the Isgur-Wise function at zero recoil [^. This effect is purely 
kinematical: it is solely due to the fact that the momentum of a quantum particle in a box is 
quantized. Because of its purely kinematical nature, we were able to correct it analytically 
and extract sensible results even though the distortions this effect gave rise to were as large 
large as 30%. For more details, please see Ref. [ES . 



2.3.4 Quenching 

As mentioned earlier, most calculations of weak matrix elements on the lattice are performed 
in the quenched approximation. This means that the fermion determinant, det(^+M), 
is set to 1 when generating the gluon configurations according to the distribution given 
in Eq. (§). Physically, this approximation corresponds to neglecting the effect of quark 
loops. The errors this approximation induces are difficult to estimate and will depend on 
the quantity studied. Experience shows, however, that the results of quenched calculations 
generally agree with experiment at a level ranging from appoximately to 20%. It should be 
noted, also, that quite a bit of progress has been made in understanding the chiral behavior 
of quenched QCD by designing and exploring a quenched version of the traditional chiral 
lagrangian. (Please see M. Golterman's lectures in the present volume for a very clear 
introduction to the subject.) 

The reason why so many lattice groups work in the quenched approximation is because, 
with present day algorithms, the CPU time required for a quenched simulation is roughly 
proportional to the number of points on the lattice while for an unquenched simulation, this 
time is roughly proportional to the number of points to the power 2.5. 

^In the quenched approximation, a virtual pion cannot be created out of the vacuum. Thus, the range 
of the strong interaction must be determined by either gluonic modes or the hghtest hadron that can be 
formed from the available valence quarks. 
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2.3.5 Matching 



On the lattice, the local vector and axial-vector currents and q''y'^j^q, where q and q' 

are quark fields, are not symmetry currents and can therefore develop a dependence on the 
cutoff through renormalization. Q To cancel this dependence and obtain physical results, 
one has to multiply the corresponding lattice results by renormalization constants which 
exactly compensate the scale-dependence of the currents. eI The renormalization constant 
for the vector current is known as Zy and that of the axial- vector current as Z^. Since 
these constants incorporate physics that lies above the cutoff, they can be calculated per- 
turbatively as long, of course, as this cutoff is much larger than Aqcd- Thanks to the 
techniques developed recently for improving the notoriously poor convergence of bare lattice 
perturbation theory [pO|] , these perturbative determinations are becoming much more reli- 



able. One is not restricted to use perturbation theory, however. One can determine Zy and 
Za non-perturbatively by tuning these constants so that hadronic matrix elements of the 
corresponding currents satisfy lattice Ward identities [^] or even some other normalization 
condition (see for example Ref. |^). For operators which are not symmetry currents in the 



continuum Ward identities are not an option and the use of a normalization condition on 
hadronic matrix elements represents a loss of predictive power. It is therefore preferable to 
normalize such operators using matrix elements of quarks and gluons. Techniques for per- 
forming such non-perturbative renormalizations of arbitrary lattice operators are currently 
being developed by G. Martinelli et al.p3|. 



In any case, whether perturbative or non-perturbative, this matching of lattice and contin- 
uum operators does introduce supplementary uncertainties whose size is determined by the 
accuracy with which the matching coefficients are calculated. 



2.4 The UKQCD Lattice 

For completeness, I now briefiy detail the parameters used for the lattice calculations whose 
results I present in the following. These calculations are performed on a X 48 lattice at 
an inverse coupling f] = 6.2 which corresponds to an inverse lattice spacing ~ 2.7 GeV 
as shown in Table |I[ We work in the quenched approximation and use the "clover" action of 
Eq. (|T9|) to describe quarks. We further rotate quark fields according to Eq. (|20| ) to obtain 
improved operators. Our leading discretization errors are thus of order O (a^a, a^). 

To keep volume errors under control, we perform our calculations with light quarks whose 
masses are in the range of the strange-quark mass (see Table |^). To obtain the corresponding 
quantities for up and down quarks, we extrapolate our results in light quark mass to the 
chiral limit (i.e. the limit of vanishing light-quark mass). Results for strange light quarks 

^For the decay B K*^, the operator whose matrix element one has to consider is aheady scale- 
dependent in the continuum because it is an effective local operator which results from integrating out the 
W. 

^These two currents do not mix with other operators on the lattice and are therefore multiplicatively 
renormalized. 
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mp ( GeV) 


0.14144 
0.14226 
0.14262 


~ 0.80 
~ 0.58 
~ 0.45 



Table 2: Kq is the light quark's hopping paramter and mp is the mass 
of the pseudoscalar meson obtained by combining this quark with its 
antiquark |jl2[. Here, rrip was used to set the scale. The errors on mp 
are dominated by the uncertainty in the scale and are on the order 10 
to 15%. 





rriQ (GeV) 


0.121 


1.90 


0.125 


1.64 


0.129 


1.36 


0.133 


1.06 



Table 3: kq is the heavy quark's hopping paramter and mq is it's 
physical mass. This mass is obtained by subtracting the 500 MeV 
of energy carried by the light degrees of freedom from the spin 
averaged mass, (3My + Mp)/4, of a vector (V) and apseudoscalar 
meson (P) obtained by combining the heavy quark with a massless 
light quark (see Ref. ||3^ for details). Here again, nip was used to set 
the scale. The errors on tuq are dominated by the uncertainty in the 
scale and are on the order 10 to 15%. 

are obtained by interpolation. 

To keep discretization errors under control, we limit our heavy quarks to have masses about 
rric (see Table ^). To obtains results which are relevant for situations where the heavy quark 
is a b, we extrapolate our results in heavy-quark mass to m^ using HQET as a guide. 

Finally, a pointlike operator such as Q{x)rq{x), where F is some combination of Dirac 
matrices, Q is a heavy quark and q a light one, usually has very poor overlap with the 
ground state of the corresponding meson since the latter has spatial extent. Thus, to improve 
the overlap it is natural to use extended operators. This usually involves "smearing" the 
heavy-quark field over a volume which is determined by a kernel K{x,x'): 

Q'ix,t)^J2K{x,x')Qi^',t), (25) 

X' 

where (^"^(x, t) is the smeared field. The results presented below were obtained using a 
gauge-invariant kernel which has an rms radius of 5.2 lattice units (for details, see ref. [p^). 
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3 Leptonic Decays ^ 



3.1 Leptonic Decays of Pseudoscalar Mesons 

In the present section, I consider the leptonic decays of a pseudoscalar meson, P, containing 
a heavy quark Q and a light antiquark q. In the Standard Model these decays are mediated 
by a boson. As mentioned in the Introduction, the coupling of the quarks to this boson is 
strongly modified by non-perturbative QCD dynamics. This modification of the quark-boson 
vertex is parametrized by a decay constant, fp, which we calculate numerically using the 
lattice formulation of QCD. 

Because the amplitude for these leptonic decays is proportional to the CKM matrix element 
VqQ {Vub for the leptonic decay of a B~ meson), one might think that a study of these decays 
could lead to a determination of VqQ. However, one can show that the branching ratio for 
these decays is suppressed by five powers of the heavy-meson's mass in the limit that this 
mass becomes very large 

BiP^iu,) (26) 

and is therefore very small. For P = and i = t~ , and for reasonable values of fs 
(~ 160 MeV) and Vub (~ 0.005), the branching ratio is on the order of 10~^. Thus it 
may be a while before these decays provide a precise determination of Vub- Nevertheless, a 
precise knowledge of the corresponding decay constants is important. These constants are 
required for describing B — B mixing llHl as well as for describing non-leptonic weak decays 
in factorization schemes ||37|| . Furthermore, a detailed study of how fp depends on the mass 
of the heavy quark, Q, provides an important tool for gauging the range of applicability of 
heavy-quark symmetry. 

To describe these leptonic decays, we must calculate the non-perturbative matrix element 

(0|W5Q(0)|P(p))=zp^/p. (27) 



As sketched in Section 2J., we obtain this matrix element from the lattice 2-point function 

E(^V75g^(x)g%gf(0)) , (28) 

X 

where the subscript s indicates that the corresponding field has been smeared as in Eq. (|25|) 
while the superscript R indicates that the field has been rotated according to Eq. (ppj). We 
evaluate this correlator for the four heavy-quark flavors described in Table ^ and for the 
three light-quark flavors described in Table H. Since the heavy-quark flavors we consider all 



^The UKQCD results presented in this section have been pubHshed in Rcf . psj . 

^°The important non-perturbative quantity here is, in fact, fs Bb where Bb is known as the "bag" 
parameter. The latter is actually currently being calculated by the UKQCD Collaboration in the limit of a 
static &-quark. 
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have masses in the range of the charm, we do not obtain Jb directly. In order to get fs, we 
must first understand how the decay constant fp varies with the mass of the heavy quark. 

HQET predicts that the quantity 

*(Mp)^(i^fM)'"°Z-/py5; (29) 

\as{mB) J ^ 

is independent of Mp up to corrections proportional to inverse powers of Mp and up to 
higher-order radiative corrections: 

HMp) ''^"^ cst + o(^,a.(Mp)) . (30) 



To evaluate <l>(Mp) we take /5o = 11 — 2/3n/ with Uf = and use for as the one-loop 
expression with Aqcd = 200 MeV. We plot our results for $(Mp) in Fig. I as a function of 
1/Mp. 

In the heavy-quark limit, <l>(Mp) is a horizontal line as Eq. ([30| ) indicates. The fact that 
$(Mp) increases quite sharply with Mp in Fig. |l| is a clear indication that corrections to the 
heavy-quark limit are quite large. To quantify these corrections we fit $(Mp) to a linear 
and a quadratic function of 1/Mp. These fits indicate that the 1/mQ-corrections to fp 
are on the order of 10-15% at the scale of the B and 27-40% at the scale of the D. This 
larger than one would expect on the naive grounds that these correction should be of order 
Aqcd/"^6 ^ 0.25/4.8 - 5% for /p and of order AqcdMc ^ 0.25/1.45 ~ 17% for fo, a 
phenomenon also observed in sumrule calculations where these corrections range from 13% 
to 22% at rriB and from 37% to 64% at mp)|]5B| . 

Having determined how fp depends on Mp, we interpolate fp in Mp to the D and extrapolate 
to the B. We present our results for the corresponding decay constants fo and fp in Table |^. 
We also provide results for fp)^ and fp^ obtained in the same way as fo and fp but with the 
light quark fixed to be the strange. We further give, for comparison, the results obtained by 
other lattice groups as well as experimental measurements when available. 

For completeness, finally, we provide in Table § the ratios of decays constants foJfD and 
fpjfp in which many systematic errors cancel. 

I should mention that there is a very vast literature on the subject of computing fp in 
the static approximation of Lattice QCD (see, for instance, the reviews of Ref.[0) and that 
a few groups are beginning to compute this same quantity using the NRQCD action of 



Eq. mm 45 



3.2 Leptonic Decays of Vector Mesons 

We now consider the leptonic decays of a vector meson, V , composed of a heavy quark, 
and a light antiquark, q. These decays are not relevant phenomenologically as a B* or D* 
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Figure 1: UKQCD data for $(Mp) plotted against inverse meson 
mass. The open symbols correspond to the decay constants obtained 
for every possible combination of our four heavy quarks mass and 
our three light antiquarks. The full symbols, on the other hand, are 
obtained by extrapolating the light-antiquark mass to zero at fixed 
heavy-quark mass. The solid line represents the linear fit to the 
chirally-extrapolated points using the three heaviest meson masses, 
whereas the dashed curve results from a quadratic fit to all four. 



meson will decay electromagnetically or through the strong interaction long before it decays 
weakly. They are interesting, however, because they provide a means of testing heavy-quark 
symmetry. HQET predicts that the quantity 
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Ref. 


Action 




/D(MeV) 


/D.(MeV) 


/d, / /d 


UKQCD|5 


1 


Clover 


6.2 


185 + ^ + 


212 + t - 1 


1.18 + 2 


APE|8[ 


Clover 


6.0 


218 ±9 


240 ±9 




BLSH 


Wilson 


6.3 


208(9) ± 35 ± 12 


230(7) ±30 ± 18 


1.11(2) ±.04 ±.02 


ELC0 


Wilson 


6.4 


210 ± 15 


227 ± 15 




ARGUS |1] 








z6 / ± Zo 




CLEOH 








344 ± 37 ± 5/ ± 42 




WA75|g 








ooo _i_ AC _i_ on _1_ AO 
161 ± 4o ± ZU ± 4o 




Ref. 


Action 




fsiMeV) 




/ /b 


UKQCD|5 


1 


Clover 


6.2 


160 + ^ + 


1 OA +6 +62 

- 5 - 9 


1.22 + \ 


APElH 


Clover 


6.0 


197 ± 18 






BLSH 


Wilson 


6.3 


187(10) ±34 ± 15 


207(9) ± 34 ± 22 


1.11(2) ± .04± .02 


ELCg 


Wilson 


6.4 


205 ± 40 




1.06 ± .04 


HEMCGCg 


W+S 


5.6 


200 ± 48 







Table 4: Heavy-light decay constants obtained with propagating, rela- 
tivistic heavy quarks. The normlization used here is the one for which 
/tt = 132 MeV. The HEMCGC calculation is unquenched and was 
performed with Wilson valence quarks and Staggered sea quarks. 



is equal to 1 up to corrections proportional to inverse powers of the mass M and higher order 
radiative corrections: 

^(M) ^'-^ \ ^o{^^,c?XM)\^ . (32) 
Here, M = (3My ± Mp)/4 and fv is the decay constant given by 

{0\q^^Q{Q)\V,e)=e^-^ , (33) 

Jv 



In Fig. 1^, our results for the quantity U are plotted as a function of 1/M. If heavy-quark 
symmetry were respected in the charm sector, we would expect all four points to lie on the 
horizontal line U = 1. The fact that U (M) is not horizontal in Fig. ^ is again an indication 
that heavy-quark symmetry is quite badly broken. To quantify this statement, we peform, as 
we did for $(Mp), a linear and a quadratic fit in 1/M to our data. We find that corrections 
to the heavy-quark limit are again on the order of 10% at the scale of the B and on the 
order of 30% at the scale of the D. 

The fact that U extrapolates close to 1 when 1/M — ^ 0-i.e. the value it is supposed to take 
in the heavy-quark limit- provides support for our description of power corrections. It also 
gives us confidence that discretization as well as other systematic errors are under control. 
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Figure 2: The quantity Lf{M) plotted against the inverse 
spin-averaged mass. Linear and quadratic fits are represented by the 
sohd and dashed curves, respectively. Also shown are the statistical 
errors of the extrapolation to the infinite mass limit. In this plot, the 
light antiquark is masseless. 



4 Semi-Leptonic Decays 

Let us now turn to the subject of semi-leptonic decays of B mesons into D or D* mesons. 
These decays are depicted in Fig. ^ As was the case for leptonic decays, the fact that the b 
and the c are confined within hadrons severely modifies their coupling to the W. 

These semi-leptonic decays are interesting theoretically. They enable one to determine the 
now famous Isgur-Wise function and to test the range of validity of heavy-quark symmetry. 
But they are also interesting phenomenologically because they can be used to measure Vet- 
The decay B D*ih' is, in fact, very well suited for this task because its rate, aside 
from being rather large, is free of (9 (AQCD/''^c)-corrections to the heavy-quark symmetry 
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prediction at maximum four-momentum transfer (Luke's theorem||46||): 



^/^J^^ duo 47r3 ' " /^'/A^ . 

where u; = vb ■ vd* is the recoil, vb{d*) are the four- velocities of the mesons and rjA are 
perturbative radiative corrections that will be described shortly. Thus, a measurement of 
the decay rate very close to the zero-recoil point yields a model-independent determination 

of 146 0] up to small, non-perturbative O ^(-^^) ^ corrections. In practice, what one does 

is to measure the differential decay rate over the whole kinematical range and extrapolate 
it to u; = 1. To be reliable, this extrapolation requires a theoretical guide, especially since 
the different experiments (ARGUS, CLEO and ALEPH) have been finding decay rates with 
rather different cu-dependences (see Figs. |T3|, |I4| and |T5|). This in turn requires a knowledge 



of the rate's dependence on u which is governed by the non-perturbative QCD dynamics 
which binds the quarks into mesons. It is here that our lattice calculation enters. 

B Diu decays are slightly less well suited for obtaining Vcb as they are helicity suppressed, 
which makes them more difficult to measure close to a; = 1. Moreover, their rate is not 
protected by Luke's theorem and should suffer larger O (Aqcd / ^c) corrections close to zero 
recoil. Neither of these problems are insurmountable, however, since the rate is far from 
negligible away from u = 1 and the O (Aqcd / ^c) corrections turn out to be parametrically 
suppressed [|]. So in the future, when better data is available, this channel should also 
provide a good means for determining Vcb- 



^The rate is down by a factor of — 1 compared to that for B D*£iy decays. 
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The matrix element required to describe B Diu and B D*£i' decays are 

and (36) 
where u = v ■ v' and e'^ is the polarization vector of the D*. 

In the limit that rric and mi, are infinite, heavy-quark symmetry reduces the six form factors of 
Eqs. (^5]) and (^) to a single universal function of the recoil, ^(cu), known as the Isgur-Wise 
function: 

h+{Lj) = hA^{uj) = hAsiuj) = hv{uj) = ^{u) 
and (37) 
h_{(jj) = h^^i^^) = . 

This symmetry further requires that the Isgur-Wise function be normalized to one at zero 
recoil 10]: 

m = 1 • (38) 

Of course, in nature, the bottom and the charm are not infinitely massive and there are 
corrections to the simple relations of Eq. (RBI) as there were to the scaling relation of Eqs. 



and (|3^). These corrections can be parametrized in terms of twelve functions (3i{uj) and 7j(c(j) 



with i = +, — , V, Ai, A2, A3, such that 

hi{u) = («, + + 7i(cu)) e(^) , (39) 

where a+ = ay = = c^As = 1 and a_ = 02 = 0. The functions Pi correspond to the 
exchange of hard gluons across the vector and axial-vector vertices. These corrections are 



perturbative and in the sequel we use the results obtained by M. Neubert in Ref. ]|47|]. In this 
reference, Neubert performs a full one-loop matching of HQET to QCD and runs the results 
at two loops. 

The second set of corrections, the 7j's, are proportional to inverse powers of the heavy-quark 
masses and rrih. They correspond to matrix elements of higher-dimension operators in 
HQET and are therefore non-perturbative and almost as difficult to evaluate as full QCD 
matrix elements. If most of these corrections cannot be neglected, heavy-quark symmetry 
obviously loses much of its predictive power. In most cases, heavy-quark symmetry is not a 
precision tool since leading power corrections are expected to be on the order of Aqcd/'"^c — 
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20%. In a few instances, however, heavy-quark symmetry is more constraining as we saw 
at the beginning of the present section when discussing the extraction of the CKM matrix 
element Vet and Luke's theorem. 

At the level of form factors, Luke's theorem guarantees that and h^-^ do not suffer 
C (Aqcd/''^c,6) corrections at zero recoil. Thus 

which means that the leading non-perturbative corrections on these form factors ought to 
be on the order of 4%-i.e., very small. 

It is important to emphasize, at this point, that to claim that the Isgur-Wise function, ^(cu), 
has been obtained from any one of the form factors, hi{u), one has to convincingly argue 
that one controls the non-perturbative power corrections, 'Ji{uj). 



4.1 The Form- Factor h+{uj) S 

The form factor dominates the rate for B Diu decays. The matrix element of Eq. (|35|) 
which h'^ parametrizes is obtained from the 3-point function 

E e-^-e-*'-y (gAscf (t/, y) c%6^(t, x) bfj,q\0, 0)) , (41) 



as sketched in Section ^Tlj (See Eqs. (psD and (^) for the meaning of subsrcript s and the 
superscript R.) 

Since the heavy-quark flavors we consider all have masses in the range of the charm, we do 
not directly obtain the form factor h'^{uj) relevant for physical B Diu decays. In order to 
get this physical form factor, we must extrapolate the mass of the initial heavy quark to rrih 
and interporlate the mass of the final heavy quark to rric. Thus, we must understand how 
h'^ depends on the mass of the initial and final heavy quarks. Since we already understand 
the dependence on mass which comes from radiative corrections, we study the quantity 

+ (42) 

where these corrections are subtracted. In this equation, equahty holds only to leading order 
in power and radiative corrections. 

In Fig. ^we plot the ratio h'^ / [l+f]'^) for four degenerate transitions. A degenerate transition 
is one in which the masses of the initial and final heavy quarks are equal. We consider these 
transitions first, because they are more constrained theoretically |3^. cl The fact that all four 



-'^^Most of the UKQCD results presented in this section will appear in Ref . [p2[ . 

^•^Because of electromagnetic current conservation, degenerate transitions have 7'''(1) = to all orders in 
l/rrib^c and h'{uj) = for all to. 
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Figure 4: The ratio h~^{uj){l + j3^{uj)) is plotted as a function of uj for 
four different degenerate transitions corresponding to four values of 
the heavy-quark mass. The data corresponding to each transition has 
it's own symbol as detailed on the plot. The different curves result 
from individually fitting each set of data to the parametrization s^^r 
(see text). The hopping parameter of the light, spectator antiquark is 
Kg = 0.14144. 



sets of points, which correspond to heavy-quark masses ranging from 1.1 GeV to 1.9 GeV, 
appear to lie very much on the same curve is a first indication that the heavy-quark-mass 
dependence of h^/{l + is very small. That all four sets of data lie on the same curve 
can actually be shown by fitting each set individually to a standard parametrization for the 
Isgur-Wise function and showing that the parameters of the resulting curves are equal 



within statistical errors. The parametrization we use is 48 



exp 



UJ 



UJ 



(43) 



where is the negative of the slope of the function ,^jv_R at = 1 and is the parameter of 
this fitting function. We acutally fit the data to s C,nr where s is a parameter which allows 



-'^"'if the heavy-quark mass dependence of (1 + (1 + /3+)) is neghgible this quantity is just the Isgur-Wise 
function ^. 



22 



us to test how well we have normalized our data. If our data is well normalized and ^nr 
is a valid parametrization for the data, then s ought to be equal to 1. For all four sets of 
data, we find s = 0.99(1) and = 1.4 where the error on s is statistical and the statistical 
errors on are the order of 0.2 but differ slighly from data set to data set. Thus, our data 
appears to be correctly normalized and independent of heavy-quark mass. 

This mass dependence can in fact be quantified by extracting the power corrections, 7"'", from 
our data. If radiative corrections are neglected, these power corrections can be parametrized 
in terms of a single universal function, g{uj), such that 

((4)^). ,44) 

where A is the energy carried by the light degrees of freedom and is approximately equal 
to 500 MeV when the light quark q is an up or down quark[^]. Furthermore, according to 
Luke's theorem, 

g{l) = . (45) 

The form factor, g{u!), can be extracted from our results for h'^ by taking ratios of the 
quantity h~^/{l + f3~^) at fixed u but different values of the initial or final heavy quark mass 
(see Eq. (^2|)). Because lattice momenta are quantized, there are, in our data, only four values 
of u for which this can be done. These four point are plotted in Fig. ^. This figure indicates 
that g{uj) is consistent with and in any case is less than one or two times 10~^ in the range 
of recoils that we can explore. Since (A/2)(l/mb + l/rric) is about 0.4 for the quarks we 
simulate power corrections to appear to be indeed very small. When combined with the 
information obtained earlier for degenerate transitions which covered a much larger range of 
recoils, this result clearly confirms that / {l + (3~^) depends very little on heavy-quark mass. 
This weak dependence on heavy-quark mass, however, can be interpreted in many ways. A 
first explanation is that power corrections to , as well as discretization errors proportional 
to the masses of the heavy quarks, are genuinely small. A second possible interpretation is 
that discretization errors are small but that the O (^A/2mfe^c) -correction happens to cancel 
against the higher order power corrections in the range of heavy-quark mass that we are 
investigating-i.e. that in this range of masses the heavy quark regime has not yet been 
reached. A third is that the power corrections and discretization errors have opposite signs 
and cancel. Since the cancellations of options two and three would have to take place over a 
rather large range of recoils and a rather large range of heavy-quark masses, they seem quite 
unlikely. So in what follows, we will adopt the interpretation that power corrections to h'^ 
are genuinely small even for heavy quark with mass close to that of the charm. 



^^These arguments are fleshed out in Ref. [Q. It must also be said that the way in which we normahze 
our data means that we measure 7'^(w) — 7^(1) and not "/^(lo). Thus, we are not sensitive to those power 
corrections which depend very weekly on lo. If these corrections happen to be large, our conclusions about 
the size of power corrections are invalid. However, because our normalization procedure subtracts some of 
the higher-order power corrections, our determination of the Isgur-Wise function as well as our determination 
of the leading O (A/2m;,^c) corrections should be all the more accurate. Again see Ref. js^ for details. 



23 



0.5 



0.0 



-0.5 



T 1 1 1 1 1 1 1 r 



O ^5 



J I I I I I I I L 



1.0 



1.1 



1.2 



Figure 5: The subleading form factor, g^u), is plotted versus u. The 
hopping parameter of the hght, spectator antiquark is Kg = 0.14144. 



The absence of heavy-quark mass dependence of the ratio /i^/(l i^l3~^) also means that to 
a good approximation we have obtained the infinite mass resuh. " Thus, we can combine 
all of our data for h^/{l + for fixed light-quark mass but all possible initial and final 
heavy-quark masses and call the resulting function an Isgur-Wise function (see Eq. (^)). 
In Fig. we plot this combined data for the case where the mass of the light antiquark is 
interpolated to the strange-quark mass. The resulting Isgur-Wise function, which we denote 
by ^s{^^), is relevant to the decays Bs — > Dgiu and Bg — > D*ih' and elastic Bs and Dg 
scattering off a photon. The different symbols on the plot correspond to different initial 
and/or final heavy-quark mass. The fact that these different symbols often lie on top of 
each other and always appear to lie on the same curve is again confirmation that the mass 
dependence of the ratio h'^ / {1 + f3'^) is very small. The solid line in Fig. ^ corresponds to a 
fit of the data to the parametrization s^nr{^), with ^jsm given in Eq. (|^). From this fit we 
find that this Isgur-Wise function has a slope of —p^ at = 1 with 

= 1.2 + 2(gtat.) + ?(syst.) (46) 

-'^^The only way /i^/(l + (3^) could not be the infinite mass results given it's apparent heavy-quark mass 
independence is if option two of the preceding paragraph were realized. 



24 




Figure 6: The ratio / (1 +/?"'") is plotted as a function of uj. The dif- 
ferent symbols correspond to different initial and/or final heavy-quark 
masses. All symbols, however, correspond to situations where the light 
antiquark has the mass of the strange. Because this ratio exhibits no 
dependence on heavy-quark mass, it is just the Isgur-Wise function 
isil^)- The solid curve depicts the result of fitting the parametrization 
s^NR to the data. 



where the sytematic error was obtained in a way described in Ref. [^. 

In Fig. 1^ we repeat the same exercise with the light-quark mass extrapolated to the chiral 
limit. The Isgur-Wise function obtained thus is the one relevant to the decays B — > Dlv 
and B D*ii> and elastic B and D scattering. This function is traditionally called "the 
Isgur-Wise function" . Performing the same fit as in Fig. |^, we find that the corresponding 
slope parameter is 

Pld = 0.9 t ^(stat.) + ^(syst.) . (47) 

This slope parameter, as well as the one of Eq. (^) is compared with the predictions of 
other authors in Table |^. Our predictions lie safely above the lower bound of Bjorken||49[ 



and below the upper bound of de Rafael and Taronpm. Our results for also agree with the 



lattice result of[51 although the details and systematics of the two calculations are different. 
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Figure 7: The ratio / [1 + f3~^) is plotted as a function of uj. The dif- 
ferent symbols correspond to different initial and/or final heavy-quark 
masses. All symbols, however, correspond to situations where the 
light antiquark is massless. Because this ratio exhibits no dependence 
on heavy-quark mass, it is just the Isgur-Wise function ^^ ^(co'). The 
solid curve depicts the result of fitting the parametrization sC,j\fji to 
the data. 



The authors of do not quote a value of for vanishing light-quark mass because of 
their poor statistics in the chiral limit. 

Also for comparison, we quote an average experimental value for the slope of the Isgur-Wise 
function compiled by Neubert ||52|| from very recent results of the ALEPH||61| and CLEOfG^ 
Collaborations as well as older data from the ARGUS Collaboration ||63i: 



Pld(e.,*.) = 0-87(12) (20) , (48) 
where the second error is theoretical and accounts for the uncertainty in the size of l/rric 



correct ions |5^. Agreement with our result is excellent (see Eq. (0))- Such good agreement, 
however, is most certainly coincidental given the size of both the experimental and lattice 
errors. 



As can be inferred from our determinations of p"^^^ and pi (Eqs. (|i7| ) and (^61) ) our results 
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Reference 






UKQCD 


0.9 + 3(stat.) + 2(syst.) 


1.2 + ^(stat.) + i(syst.) 


Bernard, Shen and Soni||51|| 




1.24(26) (stat.) (26) (syst.) 


de Rafael and Taron 


m 


p2 < 1.42 




Close and Wambach| 


a 


1.40 


1.64 


Neubert[|] 


0.66(5) 




Voloshin|p5| 


1.4(3) 




Bjorken[||] 


p2 > 0.25 


Blok and Shifman |56|| 


0.35 < p2 < 1.15 




H0gaasen and Sadzikowski[|57|| 


0.98 


1.135 


Rosner 


58] 


1.59(43) 




Burdmai] 


L 59| 


1.08(10) 




Dai, Huang and Jin|6Cf| 


1.05(20) 




Experiment (see text) 


0.87(12)(20) 





Table 5: Comparison of our lattice results for — ^^^(1) and — ^^(1) to 
the theoretical predictions of various authors and to experiment. 



suggest that decreases slightly as light-quark mass decreases. Such a decrease is consistent 
with one's intuition about the inertia of the light degrees of freedom. A very similar trend is 
also found by H. H0gaasen and M. Sadzikowski ||57||. In fact, our predictions for itself are in 
excellent agreement with theirs. Their prediction is based on an improved bag model calcu- 
lation and is an extension of earlier work by M. Sadzikowski and K. Zalewski [|64| . A similar 



decrease in slope with spectator-quark mass is observed by F. E. Close and A. Wambach||53 
though the values they quote for p^^^ and p^ are slightly larger than the ones we find. 

I wish to mention, finally, that J. Mandula and M. Ogilvie are in the process of calculating the 
Isgur-Wise function using the lattice version of HQET whose action was given in Eq. (E2|) |5| . 



4.2 The Form Factor hA^ito) 

We now turn to the subject oi B ^ D*lv decays. Here the dominant form factor is /i^i- 
This form factor is obtained from the 3-point function 

E e-^^'^e-P'-y (g%cf (t;, y) cV75&^(t, x) fef 75^^(0, 0)) , (49) 

as was from the 3-point function of Eq. (^T|). Here again, we must study the behavior 
of the form factor with heavy-quark mass to obtain the form factor relevant for physical 
B D*ii> decays since our simulation is performed with heavy-quark masses in the range 
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Figure 8: The ratio /iAi/(l + Pai) is plotted as a function of uj for 
four different degenerate transitions corresponding to four values of 
the heavy-quark mass. The data corresponding to each transition has 
it's own symbol as detailed on the plot. The different curves result 
from individually fitting each set of data to the parametrization sC,nr 
(see text). The hopping parameter of the light, spectator antiquark is 
Kg = 0.14144. 



of the charm-quark mass. We do so by considering the quantity 

-4^^^(i + 7A.M)eH, (50) 

where the equality holds only to leading order in power and radiative corrections. 

The analysis we perform of the quantity /iAi/(l + Pai) parallels the analysis of (1 + /3+) 



carried out in Section ^T^. In Fig. || we plot the ratio /iAi/(l + Pai) for four degener- 
ate transitions. Here too we fit the data for each individual degenerate transition to the 
parametrization s^nr{^), with ^nr{^^) given in Eq. (^). We find s = 1.00(2) and = 1.3 
with a statistical error on the order of 0.3 for all four data sets. This confirms the naked eye 
impression that all four data sets lie on the same curve and thus indicates that the depen- 
dence of the ratio /iai/(1 + Pai) on heavy-quark mass is very small in the range of masses 
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we are considering. Moreover, the fact that the curves on which all of these points lie are 
almost the same as those found for + (3^) is a first indication that the spin component 

of the heavy-quark symmetry is unbroken in this case even in the region of the charm-quark 
mass. 

As we did for + /S"*"), we can try to quantify the heavy-quark mass dependence of the 

ratio + /?yii) by measuring the power corrections to this ratio. If radiative correc- 

tions are neglected, the power corrections to can be parametrized in terms of the three 
univsersal functions g{uj), rji^uj) and g*{uj): 



where g{io) is the same form factor that appears in the power corrections to in Eq. 
Here again, Luke's theorem requires that the O {K/2mi,^^ power corrections to Ha^ to vanish 
at = 1. This means that g*{uj) must vanish at = 1: 

(7*(1)=0. (52) 

?7, however, is unconstrained because it does not appear in the expression for power correc- 
tions aX uo = 1. 

We extract the functions which multiply the expansion parameters A/(2mb^c) in Eq- (|5l|) by 
taking ratios of the quantity + (3ax) at fixed uj and light-quark mass but different 

values of the initial or final heavy-quark mass. We plot our results in Fig. ^ Both these 
functions are consistent with and are less than one or two times 10~^ in the range of recoils 
that we can explore. Since these functions are themselves multiplied by the small expansion 
parameters A/(2mc) and A/{2mb) which are on the order of 0.2 for the quarks we simulate, 
the power corrections to Ha^ are very small (on the order of 4% or less) in the explored range 
of recoils. Together with the information provided by the fits in Fig. ^, this results indicates 
that power corrections to Ha^ are very small for all u, up to the caveats mentioned in Section 
4?T| in the discussion on power corrections to h^. 



The fact that hA^/{l + PaJ depends very little on heavy-quark mass means, as it did for the 
ratio /i^/(l + 13~^), that our results are to a very good approximation infinite mass results. 
So again we can combine all of our data for a given light-quark mass to obtain the Isgur- 
Wise functions C,s and C,u,d which are in principle the same functions as the ones given by 
the ratio h~^/{l + We plot these form factors as functions of the four- velocity recoil in 



Figs. M and 11. The solid curve in both these plots result from fitting our data to the 



parametrization s^nr- These fits, give the following slope parameters: 

pI = 1.2 + ^(stat.) + ^(syst.) (53) 



^^Here too, the manner in which we normahze our data means that we measure 7^j(a;) — 7^j(l) and not 
7^j(aj). The consequence of this are the same as those mentioned regarding power corrections to h'^. 
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Figure 9: Power corrections to the form factor h^^. Plotted are 
the functions which multiply the expansion parameters A/(2mb^c) in 
Eq. (pT]). The hopping parameter of the light, spectator antiquark is 
Kg = 0.14144. 



and 

pIa = 0-9 t 5(stat.) t ?(syst.) . (54) 



These slopes are precisely the ones we found for the ratio h^/{l + (see Eqs. ( ^^ and 
(|47|)). The Isgur-Wise functions obtained from and are in fact equal over the whole 
range of recoils to within less than 4%, as can be seen in Fig. ^ where we plot the ratio 

hAA^ + P+)/h+{l + (3A,). 

What has emerged is a very consistent picture of i? — >■ Diu and B D*(!.v decays in which 
heavy-quark symmetry is surprisingly well satisfied even though the heavy quarks with which 
we work have masses in the range of the charm-quark mass. This is in stark contrast with the 
results for the decay constant fp presented in Section |3.1| where we found that corrections to 



the heavy-quark limit were on the order of 30% for these same heavy quarks. What seems to 
be happening here is that the protection from O (^A/2mb c) -corrections that Luke's theorem 
provides at zero recoil appears to extend over the full range of recoils so that corrections 
which one would naively expect to be on the order of (A/2)(l/mc + l/m?,) ~ 30 — 40% for 
the quarks we are studying turn out to be on the order of a few percent. 
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Figure 10: The ratio + PaJ is plotted as a function of u. 

The different symbols correspond to different initial and/or final 
heavy-quark masses. All symbols, however, correspond to situations 
where the light antiquark has the mass of the strange. Because this 
ratio exhibits no dependence on heavy-quark mass, it is just the Is- 
gur-Wise function C,si^^)- The solid curve depicts the result of fitting 
the parametrization s^nr to the data. 



4.3 Extraction of Vcb 

When power corrections and radiative corrections for u; > 1 are neglected, the differential 
decay rate for B —>■ D*lv decays is given by 



1 dV{B D*iu) ^ Gl\Vc 12 
- 1 duj ~ Att^ 



{niB - mD*ym%,r]\{l + 5i/m2) ^u,di^) ^ (55) 



where Si/m^ stands for the power corrections to h^^ at c<j = 1 which have been the object 
of much controversy lately [|^, Having determined ^u,d with our lattice calculation, the 
only unknown left in Eq. (|55D is the CKM matrix element | Vet \ ■ Thus, a fit of the theoretical 
expression of Eq. (pSj) to an experimental measurement of the rate immdediately yields a 



measurement of \Vcb\- This is what we do in Fig. 13 where we use very recent data obtained 



by the CLEO collaboration p2|. The value of \Vcb\ that this fit gives is 
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Figure 11: The ratio + PaJ is plotted as a function of u. 

The different symbols correspond to different initial and/or final 
heavy-quark masses. All symbols, however, correspond to situations 
where the light antiquark is massless. Because this ratio exhibits no 
dependence on heavy-quark mass, it is just the Isgur-Wise function 
^u,d{i^)- The solid curve depicts the result of fitting the parametriza- 
tion sC,NR to the data. 



IK 



cb\ 



0.037 



+ 1 

- 1 



2+4 
2 - 1 



0.99 



r]A J l + 5i 



(56) 



where the first set of errors is due to the experimental uncertainties, the second set of errors 
to the statistical errors in our determination of the Isgur-Wise function and the third to our 
systematic errors. A similar fit to ALEPH data||61| gives (see Fig. |1^ 

: + 2 + 4 /l.53ps /0.99\ 1 



ch\ 



0.042 



- 2 



l/m2 



(57) 



and to ARGUS data||63| (see Fig. 
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Figure 12: The ratio of the Isgur-Wise functions obtained from 
and /zai is plotted as a function of u. Only degenerate transitions are 
considered. The different symbols correspond to different values of the 
heavy quark's mass. All symbols, however, correspond to situations 
where the light antiquark has a hopping parameter Kg = 0.14144. 
The points are very close to one, as they should be if heavy-quark 
symmetry is respected. 



These results for \ Vcb\ must not be taken too literally because the experimental measurements 
are binned according to slightly biased estimators of the recoil and because we have neglected 
small l/rrib^c and radiative corrections for > 1. Neverthless, it is clear from Fig. O, Fig. P 



and Fig. |T5| that our prediction is consistent with experiment and favors slightly the data of 
CLEO and ALEPH over that of ARGUS. 



4.4 The Form Factor hv{io) 

We now briefly turn to the form factor hy defined in Eq. (|37|). This form factor is interesting 
because it is not protected by Luke's theorem as are and Hai- This means that one would 
expect power corrections to this form factor to be more in line with naive expectations. 

We have very prelimary results for the form factor hy and it does indeed seem to exhibit 
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Figure 13: x^-St of the expression + ^i/m^) to the 

experimental data for this quantity obtained by the CLEO Collab- 
oration In this fit, the function ^^,^(0;) is fixed to our lattice 



prediction, i.e. to with given by Eq. (|47|). The fit parameter 
is \Vcb\ TjA (1 + ^i/m^)- The UKQCD data and curve are just those of 
Fig. 1^ appropriately rescaled by the fit parameter. The experimental 
data assumes tbo = 1.53(9)ps and tb+ = 1.68(12)ps. 



power corrections on the order of 20 to 40% depending on the heavy-quarks considered. It 
is unclear, however, what fraction of these corrections are true l/m;, c-corrections and what 
fraction are amq discretization errors. The problem here is that we cannot subtract these 



discretization errors as we did for and /i^i (see Ref. for details), because there is no 
normalization condition for hy^ 

So instead of presenting results which may suffer from large discretization errors, we will 
present a framework in which the form factor hy may be analyzed once these errors are 
controlled. As we did for these two other form factors, we should study hy with the radiative 



^^For the case of /i+ and Hai , Luke's theorem guaranteed the these form factors would be 1 at a; = 1 up 
to calculable perturbative correctfons and small non-perturbative corrections O I {K/rrib^c 
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Figure 14: Same fit as in Fig. [T^ but to data of ALEPH[^ wliicli 
assumes tbo = 1.53(9)ps. 



corrections taken out: 



hv{uj) 



(59) 



wliere 7y labels power corrections. Tliese corrections can be parametrized by tlie tliree 
universal functions g, g* and r] that we encountered in Section [4.1| and Section [4.2| . In the 
absence of radiative corrections, 



7v 



c,b, 



(60) 



Then, by considering ratios of the quantity /iy/(l + Pv) for fixed u and light-quark mass, 
but different heavy-quark masses, we would extract the functions g{uj) -|- (1 — 2r]{uj)) and 
g*{uj) + 1. We would expect both these functions to be of O (1) since there is no symmetry, 
here, which forbids the appearance of power corrections. It is interesting to note that this 
expectation is consistent with the statement that g*{uj) must vanish at a; = 1 (see Eq. (|5^) ). 

Having obtained these functions, we would combine them with the results for the power 
corrections to and h^^ of Figs. ^ and ^ and solve the resulting system of equations at 
fixed uj for the subleading universal functions g, g* and rj. 
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Figure 15: Same fit as in Fig. |T3| but to data of tlie ARGUS 
Collaboration ||63| which assumes tbo = 1.53(9)ps. 



We have actually carried out this whole procedure on our preliminary data and find results for 
g, g* and i] which are very much in line with what one can infer from sumrule calculations H. 



5 Penguins on the Lattice 
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We now turn to the study of -B ^ K*'y decays. These decays occur through the quark level 
process b ^ s'y depicted in Fig. |l^ which mediated by a flavor-changing-neutral-current 



(FCNC). Such FCNC processes are interesting because they are fobidden at tree level in the 
Standard Model and only occur at one or higher-loop order. Their study therefore provides a 
means for testing the details of the Standard Model. And because even at lowest order they 
are sensitive to the presence of new particles that may appear in the loops of their diagrams, 
they may even give us a handle on physics beyond the Standard Model at comparatively 
low energies. In fact, bounds on the 6 — 57 branching ratio have already been used to place 
constraints on supersymmetric as well as non-supersymmetric extensions to the Standard 



^^Much the UKQCD data presented in this section have appeared in Rcf. 
version. 



or will appear in its revised 
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u, c, t 



Figure 16: Example of a penguin diagram that contributes to the 
decay b ^ s'j. 



Model. A comprehensive review of these results can be found in Ref . \p7\\ . 

Another interesting aspect of 6 — 57 decays is that they may permit a measurement of the 
poorly determined CKM matrix element Vts- 

Experimentally, b ^ s'j decays are difficult for the same reasons that make them interesting 
theoretically. They occur at one loop and their rate is suppressed by two powers of the Fermi 
constant Gp: they are rare decays. Nevertheless, the CLEO Collaboration has been able to 
measure the branching ratio for the exclusive channel B — > K*'y. It found that this decay 
has a branching ratio 

Bn {B K*-f) = (4.5 ± 1.5 ± 0.9) x 10"^ . (61) 

In order to compare this very pretty experimental result with the predictions of the Standard 
Model, one has to compute the long-distance contributions of the strong interaction which 
are, of course, non-perturbative. These long-distance contributions are given by the hadronic 
matrix element {K* \ S(T^i,g'^(l + 75)6 \B) which can be parametrized by three form factors [p^ 

{K*{k,e)\-sa,^q^^^^±^b\B{p))=j:c;T,{q% (62) 
^ i=i 

where k and e are the momentum and polarization vectors of the K*; p is the momentum of 
the B] q = p — k; and 

Cl = 2£^,,,eVfc^ (63) 
Cl = ie^{ml - m\.) - it ■ q{p + k)^, (64) 

Cl = ^e-gL--^^^(p + /c),). (65) 
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On the lattice, these form factors are obtained from a 3-point function in very much the 
same way as were the semi-leptonic form factors of Section ^. The viabihty of this particular 
calculation was first demonstrated by Bernard, Hsieh and Soni in Ref. . 

As the photon emitted is on-shell, the form factors need only be evaluated at g^=0. In this 
limit, 

T2(g2=0) = Ti(g2=0), (66) 

and the coefficient of T^{q'^=0) vanishes. Hence, the branching ratio can be expressed in 
terms of a single number which we take to be Ti{q'^=0). 

As in our studies of leptonic and semi-leptonic decays of B mesons, the simulation is per- 
formed with the four heavy quarks listed in Table 0. Since these quarks have masses around 
that of the charm, results have to be extrapolated to rrib. Here, however, the dependence of 
the relevant form factors on heavy-quark mass is not as straighforward as it was for leptonic 
decay constants and semi-leptonic form factors because the light degrees of freedom in this 
decay have momenta comparable to the mass of the 6-quark in large sections of phase space. 
Nevertheless, in a region around = (m^ — rriK'Y heavy-quark symmetry will apply. In 
that region, one finds that Ti(g^) and T^i^q^) scale according to[[fO[|: 



and 



Ti (q^ ^ qiax; ^kA = aoV^ x [a,(mp)]-2M' [ i + + ^ + . . . ] 

\ TTlp Tfl p J 

T2 (q' ^ gL.; mp; m^,) = bo^= x [^.(mp)]-^//^" [ 1 + + _^ + . . . ] , 
^ ^ y^'mp \ mp nip J 



(67) 



(6^ 



with the same one-loop running coupling ag and the same [3q as for the leptonic case (see 
comments after Eq. (|30|) ). 

The problem now is to find a way to use these scaling relations to get Ti (g^=0;mp; rriK*)- 
The first approach is to: 

1. extrapolate the direct measurements of T2 (g^ = q'^axi ^k*) to mp = niB] 

2. use the results for this form factor at fixed mp ~ m^) to make an educated guess about 
its behavior in at m^; 

3. use this behavior to extrapolate the value of T2 (g^ = q^axi ^b] ttik*) obtained in step 
1 to g2 = 0. 



This method has the advantage that it only requires knowledge of T2's q^ dependence at m^. 
Its disadvantage, however, is that one has to extrapolate T2 over a wide range of momentum 
transfers from q^ = q'^^x to = 0. Since one does not exactly know the functional behavior 

^°This approach can also, of course, be followed with Ti instead of T2. 



38 



0.6 



T2(q^), pole fit 



0.4 




0.2 



/Ci=0. 14144 
/Cj^=0. 13300 



K =0.14226 



0.0 



-0.2 



-0.1 0.0 



0.1 



Figure 17: T2 plotted as a function of g^. Even though T2 is consistent 
with pole dominance (solid line), it is clear that T2 could also be 
constant. The dotted lines represent the 68% confidence levels of the 
fit to the pole form. 



of T2 in g^, this extrapolation leads to rather large uncertainties. Bernard, Hsieh and Soni 
make a rather good case for the use of pole dominance ||69||, but it is difficult to exclude the 
possibility that T2 might be a constant on the basis of the data alone (see Fig. O taken from 
Ref.g). 

Step 1 is performed by extrapolating the quantity 



f2{mp) = T2(g^„^).p^ f oi,(jnp) \ 



2//3o 

(69) 



linearly and quadratically in 1/mp to l/rriB (see Eq. (|68D). This extrapolation is shown in 
Fig. ||. It yields 

T2{q'=ql,,;mB;mK^) = 0.269^^^ ± 0.011 , (70) 

taking the quadratic fit as the best estimate and the difference between the central values 
of the linear and quadratic fits as an estimate of the sytematic error (the second error in 
Eq. m). 
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Figure 18: T2(mp) is extrapolated lineary (dotted line) and quadrat- 
ically (solid curve) to mp = niB- 



Now, if pole behavior is assumed for T2(g^; m^; m^^.), this result for T2(g^=g^^^; m^; m^^*) 
implies (see Eq. (|66D) 

T,{q^=0;mB;mK*) =Tt\q^=0-,mB-,mK*) = 0.112tlt\l (71) 

where the first error is statistical and the second is the systematic error obtained by combin- 
ing the variation of the pole mass within its bounds and the systematic error from Eq. (ffOl). 
The particle exchanged here is the 1^, Bgi state whose mass has not yet been measured. The 
authors of Ref.[p^ are nevertheless able to estimate this mass to be m^^j = 5.74 ±0.21 GeV 
by using HQET. 

If one assumes, on the other hand, that T2(g^; m^; rrix*) is constant, then Ti(g^=0; m^; rrix*) 
is immediately given by Eq. ([701). 



The second approach to obtaining Ti(g^=0; m^; rriK*) is to translate the scaling relations of 
Eqs. (|67D and (|68D to scaling relations for Ti^2 (i?^ = 0; mp; rriK*)- This requires an assumption 
about the dependence of the form factors for all mp. As was pointed out by As. Abada 
in Ref.[7l|, it is inconsistent to assume that both form factors behave according to a pole 
dominance form. Indeed, if one combines this assumption with the scaling relations of 
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Eqs. (0) and (^), one finds that the form factors scale at g^=0 as Ti(0) ~ mp^"^ and 
r2(0) ~ vtip^'^ in clear contradiction with the fact that T'i(O) = T2(0) (Eq. (|66D). Thus, only 
one of the form factors can obey pole dominance. If one assumes that Ti(g^) does, which is 
consistent with T2(g^) being constant, then Eq. (|67D implies that Ti(0;mp) scales according 
to: 

Ti (g2=0; mp- mx*) = c^rup^''' x [a,(mp)]-2/^« [ i + ^ + ^ + . . . ) . (72) 
^ ' \ mp mp ) 

If, on the other hand, pole dominance for T2{q^) is assumed, which is consistent with Ti(g^) 
having a dipole form, then by Eqs. (|68|) and ( |66D Ti(0; mp) scales as: 

Ti (g2=0; mp; m^.) = domp^/' x [a,(mp)]-2/'3« ("i + A + + . . . (73) 

To obtain these scaling relations one has to expand the pole masses and q^ax inverse powers 
of mp which is reasonable in the limit of large mp. The advantage of this second approach 
is that it does not require an extrapolation over a wide range of momentum transfers. Its 
disadvantage is that it requires one to assume that the form factors has a certain form for 
all mp in the range from ~ mpi to mp. 

To implement the scaling relation of Eq. (|7^), we construct the quantity 

/ / / \ \ 2//3o 



y mp \as{mB) J 

^ — 1/2 

as we did for the scaling of T2. Extrapolating T[" (mp) quadratically in 1/mp to 1/mp 
then gives (see Fig. |^) 

Ti (g2=0; mp; m^*) = 0.159 t It . (75) 

There is a problem here, however. One can also obtain the scaling relation of Eq. (|7^) by 
combining Eqs. (|66|) and (|6^) with the assumption that T2(g^,mp) is a constant function 
of for all mp in the range of mp, to mp. But if one does this, one has to conclude that 

^ ^ —1/2 1— I 

T2{mp) = (mp) in contradiction with the results of Figs. |19| and |1^. So it appears that 

the preliminary data of the UKQCD Collaboration is inconsistent with the scaling relation 
of Eq. (IT2D and therefore with the assmptions that led to it. 

The extrapolation of Ti (g^=0; mp; m^*) according to Eq. ([7^) is performed, on the other 
hand, by extrapolating the quantity 

fr'"(mp) = T,(f=0) f^f' ( (76) 

Vmp/ \as{mB) J 

quadatrically in 1/mp as shown in Fig. This extrapolation gives 

Ti (g2=0; mp; m^*) = 0.123 t H ■ (77) 
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Figure 19: The solid curve represents the quadratic extrapolation of 

^ — 1 /2 

{mp) to rrtx* /"^s. The dotted curves correspond to the statis- 
tical errors on the fit. 



method 


Ti {q^=0;mB;mK*) 


T2{qliax) ~ rri'^/'^, T2{q^]niB) constant 


0.2691^'^ ± 0.024 


Uqi,,)^m-y\T2{q'-mB) pole 


0.1121^ ±0.018 


Ti(0) ~ m^^"^ (T2(g^;mp) constant) 


0.159 li]3± 0.015 


TM-rnf'^ {T,{q';mp) pole) 


0.123 ±0.010 



Table 6: Preliminary UKQCD results for Ti (g^=0; m^; m/^.) as ob- 
tained from the various assumptions described in the text. 



In Table I summarize the different values for Ti (g^=0; m^; m^-*) obtained above. I add 
to these results a systematic error of asamq ~ 8%, where rriQ is the bare mass of our 
heaviest heavy-quark, to account for possible discretization errors as discussed in the section 
on Symanzik improvement. The consistency of the T2-pole results plus the fact that the 
scaling relation of Eq. (^) appears to be inconsistent with the preliminary lattice results of 
the UKQCD Collaboration together seem to indicate that a value for Ti{0,mB) of around 
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Figure 20: The solid curve represents the quadratic extrapolation of 
T{" (mp) to mxt/mB- The dotted curves correspond to the statis- 
tical errors on the fit. 



0.11-0.12 is favored. However, since these results are still preliminary and do not belong to 
me, I shall not venture a final number. Instead, in Table |^ I provide the results obtained by 
other lattice groups for comparison. 

Finally, for comparison with experiment, it is useful to convert Ti(0) into a value for the 
dimensionless hadronization ratio 



R 



K* 



\mb/ \ 



m 



K* 



m. 



\Ti{q'=0)\' 



(78) 



Many of the theoretical uncertainties which arise in relating the amplitudes for these decays 
to their branching ratios cancel in this ratio. The values for Rx* obtained from the results 
for Ti(0) summarized in Table |^ are given in Table |^. These values for Rx* are independent 
of whether one considers — K*^'y or — >■ K*^'y decays. I have chosen rub = 4.87 GeV 
for consistency with Ref.|7^ and used the 1994 Particle Data Book[^ for all other masses. 
Also given in Table 13 is the experimental result. This number is obtained by taking the ratio 
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Ref. 


Action 


/5 


method 


APE| 


72 1 


Clover 


6.0 


Ti(0) ~mp^/^ 


Ti(0) ~mp^/^ 




0.21(4)(2) 


0.12(2)(3) 


BHS[3^| 


Wilson 


6.3 




Tiiq^-.niB) pole 


0.10 ±.01 ±.03 


LANL 


74 1 


Wilson 


6.0 


Ti(0) ~mp^/^ 


Ti(0) ~mp^/^ 




0.25(2) 


0.09(1) 



Table 7: Results for Ti{q'^=0;mB;mK*) obtained by other lattice 
groups making use of a variety of methods described in the text and 
summarized in Table ^ All of these results are preliminary except for 
those of BHS. The LANL group also gets, using the method of BHS, 
Ti(g2=0;mB;m^,) = 0.10(1). 



method 


Rk* 


T2{qmax) ~ m'^/'^, T2(g^;ms) constant 


(34 1 ^±3 ±5)% 


T2{ql,ax) ~ T2(g2;mp) pole 


(5.9 ^°:^±1.9)% 


Ti(0) ~ nip^'^ {T2{q^;mp) constant) 


(12 t 5 ±2)% 


Ti(0) ~mp=^/' (T2(g2;mp) pole) 


(7.1 1 ^:^±i)% 


CLEO (see text) 


(19 ± 13)% 



Table 8: Comparison with experiment of preliminary UKQCD results 
for Rk* obtained from the results for Ti (g^=0; m^; rriK*) of Table ^. 



of CLEO's measurement of the branching ratio BTZ {B K*'~f) given in Eq. dll) to their 
measurement of the inclusive branching ratio BTZ (b 57) = 2.32 ± 0.51 ± 0.29 ± 0.32|7^ . 
Both experimental and lattice uncertainties are so large at this early stage that it is difficult 
to draw any firm conclusion from a comparison. 

Before ending this discussion of radiative B decays, I would like to make a brief comment 
on corrections to the heavy quark limit. Fig. ^ shows that T2 suffers 13% corrections at the 
scale of the B and 27% corrections at m^. This is very much in line with what we found for 
leptonic decays. 



6 Conclusion 



Lattice QCD studies are already providing information about the strong interaction effects 
in the weak decays of S-mesons which is of fundamental phenomenological and theoretical 
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importance. 

As far as phenomenology is concerned, we have seen that lattice studies of leptonic decays of 
S-mesons have lead to predictions for the decay constant fs required for describing B — B 
mixing as well as non-leptonic decays in factorization approximations [p7|] . These predictions 
are well summarized by the statement 

/b = 180 ± 40 MeV . (79) 

We have also seen that lattice simulations can be used to determine the form factors required 
for guiding the extraction of the CKM parameter \Vcb\ from experimental measurements of 
the differential decay rate for B D*iv decays. In Ref. [^] these form factors are further 
used to predict various semi-leptonic i?-meson decay rates. Finally, we saw that the lattice is 
beginning to make predictions for the form factors relevant to the rare decay B K*'~f and 
for the corresponding hadronization ratio Rk* ■ Because rare decays are sensitive, low-energy 
probes for physics beyond the standard model, these predictions are very important. 

On the theoretical side, the fact that quark masses are adjustable in lattice calculations 
enables one to trace out precisely the dependence of various quantities on heavy-quark mass 
and hence probe the range of applicability of heavy-quark symmetry and test the validity 
of HQET. When studying leptonic decay constants we found that power corrections to the 
heavy-quark limit were on the order of 10-15% at the scale of the 6-quark and 30-40% at rric. 
This is significantly larger than one would expect on the grounds that these corrections are 
proportional to Aqcd/"'^;) — 5% and Aqcd/'"^c — 17%, respectively, but is consistent with 
other theoretical determinations! 



When studying semi-leptonic B Diu and B — > D*ii' decays, on the other hand, we found 
that the form factors h'^ and h^-^ suffer power corrections which are much smaller than one 
would expect on those same grounds. We took this to indicate that the protection from 
C (AQCD/''Tic,6)-corrections at zero recoil that Luke's theorem provides appears to extend 
over the full range of recoils. It also meant that two independent determinations of the 
Isgur-Wise function could be obtained from our results for /i"*" and /i^j. The two Isgur- 
Wise functions found in this way were identical indicating that the spin component of heavy 
quark symmetry is nearly unbroken in this particular situation. This procedure for obtaining 
the Isgur-Wise function was repeated for two values of the light-quark mass: rriq = and 
viq = rris. The slopes or the corresponding Isgur-Wise functions, ^u,d and ^s, are at u = 1: 

e:,dU=i = - [0.9 t ^(stat.) t ^(syst.)] (80) 

and 

e:U=i = - [1.2 + ^(stat.) + lisjst.)] . (81) 

We were also able to extract some of the form factors which appear at order l/m^ c in the 
description of the semi-leptonic decays, thereby probing some of the more intricate details 
of HQET. 
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When studying the radiative decay B — > K*'j we found that T2 defined in Eq. ( |69D suff'ered 
power corrections on the order of 10% at mf, and 30% at m^, very much in fine with leptonic 
decays. This information, however, does not tell us how the process at q^=0 scales with 
the mass of the initial heavy quark. To determine this scaling one has to understand the 
dependence of at least one of the relevant form factors as discussed in Section ^ Since 
the precision of the lattice prediction is vitally dependent on understanding this scaling, it 
is important to settle this issue soon. 

It is also important, in the near future, to improve calculations of the S-parameter for B — B 
mixing and to perform systematic studies of semi-leptonic B p{TT)iu. These processes 
can be approached with the methods described above and are important for determining 
the CKM matrix elements \Vtd\, \Vts\ and \Vub\, respectively. To permit an alternative de- 
termination of iV^bl, the UKQCD Collaboration is furthermore undertaking a study of the 
semi-leptonic decays of baryons, such as the Af,, which contain one heavy quark and two light 
quarks. In the process, many more tests of heavy-quark symmetry will be performed. Some 
time and effort should also be spent on non-leptonic decays. B J/'^Kg or B 7t~^7t~ 
decays, for instance, are important for understanding CP violation in the bottom quark 
sector. The stumbling blocks here are both technical and conceptual. Technically, present 
day lattices are just too small to separate the hadrons in the final state and conceptually, 
how to determine the relevant final-state phase shifts is not yet understood |76[ . 



The results presented here are, for the most part, products of a first generation of calculations. 
The errors quoted on these quantities will therefore decrease in the months and years to come 
as we control better the various sources of systematic uncertainties and design more efficient 
lattice actions and faster algorithms. When working with heavy quarks, aside from un- 
quenching of course, the challenge will be to reduce discretization errors. This will involve 
extrapolating results to the continuum limit by performing high statistics simulations on 
lattices with different of lattice spacings. It will also involve designing actions, such as 
the Clover or even perfect actions mentioned above, which converge more rapidly to the 
continuum limit. One further expects effective theories, such as NRQCD or the lattice 
variant of HQET, to play a growing role in the description of heavy-quark decays. In any 
event. Lattice QCD has already shown itself to be a valuable quantitative tool for strong 
interaction physics and it will undoubtedly provide increasingly precise and varied results in 
the future. 
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